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The QCD critical point can be found in heavy ion collision experiments via the non-monotonic 
behavior of many fluctuation observables as a function of the collision energy. The event-by-event 
fluctuations of various particle multiplicities are enhanced in those collisions that freeze out near the 
critical point. Higher, non-Gaussian, moments of the event-by-event distributions of such observ- 
ables are particularly sensitive to critical fluctuations, since their magnitude depends on the critical 
correlation length to a high power. We present quantitative estimates of the contribution of critical 
fluctuations to the third and fourth moments of the pion, proton and net proton multiplicities, as 
well as estimates of various measures of pion-proton correlations, all as a function of the same five 
non-universal parameters, one of which is the correlation length that parametrizes proximity to the 
critical point. We show how to use nontrivial but parameter independent ratios among these more 
than a dozen fluctuation observables to discover the critical point. We also construct ratios that, if 
the critical point is found, can be used to overconstrain the values of the non-universal parameters. 



I. INTRODUCTION AND ILLUSTRATIVE 
RESULTS 

One of the main goals of heavy ion collision experi- 
ments is to map the phase diagram of QCD. The second- 
order critical point at which the first-order transition 
between hadron matter and quark-gluon plasma (QGP) 
ends is one of the distinctive features of the phase di- 
agram. We currently do not have a systematic way of 
locating this point from first principles as model and lat- 
tice calculations face many challenges and much work 
still needs to be done in order to overcome them. (For 
reviews, see Refs. [iHS].) In the meantime, if the critical 
point is located in a region accessible to heavy-ion col- 
lision experiments, it can be discovered experimentally. 
Experiments with this goal are underway and planned 
at the Relativistic Heavy Ion Collider (RHIC) at the 
Brookhaven National Laboratory (BNL) and at the Su- 
per Proton Synchrotron (SPS) at CERN in Geneva [T(J- 
[T5] . It is therefore important to define, evaluate the util- 
ity of, and select experimental observables that will allow 
us to locate the critical point, if it is located in an exper- 
imentally accessible region. 

In heavy-ion collision experiments, the center of mass 
energy ^/s is varied, thus changing the temperature and 
chemical potential of the produced matter and in this 
way scanning the phase diagram. The observed collec- 
tive flow of the produced matter at RHIC strongly sug- 
gests the production of a strongly-coupled quark-gluon 
plasma [Hj. As the QGP expands and cools, it follows a 
path on the phase diagram that is characterized by ap- 
proximately constant entropy density to baryon number 
density ratio until freeze-out, after when there are no fur- 
ther interactions that change the multiplicities of hadron 
species. When these particles are then detected, they 
give us information about the state of the matter at the 
freeze-out point. Therefore, in order to see the effects of 



the critical point on observables, one should try to get 
the freeze-out point as close to the critical point as pos- 
sible by varying the collision center of mass energy, ^/s. 
Decreasing ^/s decreases the entropy to baryon number 
ratio, and therefore corresponds to increasing the baryon 
chemical potential /is at freeze-out. 

Present lattice calculations evade the fermion sign 
problem in different ways that all rely upon the small- 
ness of ^s/(3T). Although each is currently limited by 
systematic effects, and they do not give consistent guid- 
ance as to the location of the critical point, all present 
lattice calculations agree that it is not found at hb < T, 
where the calculations are most reliable p!5HT5] . For this 
reason, experimental searches focus on collisions which 
freezeout with /is > 150 MeV. The upper extent of the 
experimentally accessible region of the phase diagram is 
determined by the largest freezeout /is at which colli- 
sions still have a high enough -^/s that the matter they 
produce reaches temperatures in the transition region. 

Upon scanning in ^/s and thus in /i^, one should then 
be able to locate (or rule out the presence of) the critical 
point by using observables that are sensitive to the prox- 
imity of the freeze-out point to the critical point [HI |^ . 
For example, for particles like pions and protons that 
interact with the critical mode, the fluctuations in the 
number of particles in a given acceptance window will 
increase near the critical point as the critical mode be- 
comes massless and develops large long- wavelength corre- 
lations. As we vary -^s, therefore, if the freeze-out point 
approaches the critical point, we would see an increase in 
the fluctuations in the number of those particles which 
interact with the critical mode. These fluctuations would 
then decrease as we move away from the critical point. 
(This is true for any observables which are sensitive to 
the proximity of the critical point to the point where 
freeze-out occurs.) Hence, a characteristic signature of 
the critical point is the non-monotonic behavior of such 
variables, as a function of -y/s [111 HH]- Another way to 
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change the freeze-out point is by changing the size of the 
system by varying the centrahty of the coUisions, since 
larger systems freeze-out later and hence at somewhat 
smaller temperatures. 

In this paper we describe how to use the increase in 
fluctuations of particle numbers near the critical point as 
a probe to determine its location. The way one character- 
izes the fluctuations of an observable is by measuring it in 
each event in an ensemble of many events, and then mea- 
suring the variance and higher, non-Gaussian, moments 
of the event-by-event distribution of the observable. The 
contribution of the critical fluctuations to these moments 
is proportional to some positive power of ^, the correla- 
tion length which, in the idealized thermodynamic limit, 
diverges at the critical point. In reality, ^ reaches a max- 
imum value at the critical point but does not diverge 
because as it cools the system spends only a finite time 
in the vicinity of the critical point. The system also has 
only a finite size, but it turns out that the finite time is 
a more stringent limitation on the growth of the correla- 
tion length [2ni m] ■ Estimates of the rate of growth of ^ 
as the collision cools past the critical point (which take 
into account the phenomenon of critical slowing down) 
suggest that the maximal value of ^ that can be reached 
is around 1.5 — 3 fm [21 23 , compared to the natural 
~ 0.5 fm away from the critical point. Higher moments 
depend on higher powers of ^, making them more favor- 
able in searching for the critical point [21] . In this paper 
we consider the second, third and fourth cumulants of 
particle multiplicity distributions for pions and protons. 
We also consider mixed pion-proton cumulants, again up 
to fourth order. 

Our goal in this Introduction is to provide an illus- 
trative example of one possible experimental outcome. 
In Section LA we define the observables that must be 
measured at each ^/s. In Section LB we suppose that 
the critical point is located at /is = 400 MeV and then 
guess how the correlation length ^ at freezeout will vary 
with the chemical potential /is, and hence with y^, in 
a heavy ion collision program in which the beam energy 
is scanned. In Section I.C we plot results for how seven 
of the observables that we define will vary with ^b, if 
the guess for £,{1^3) that we have made for illustrative 
purposes were to prove correct. In Section II we provide 
the calculation of all the observables that we define, as a 
function of ^, the proton and pion number densities, and 
four nonuniversal parameters that must ultimately be ob- 
tained from data. In Section III we construct ratios of 
observables that allow us to measure four combinations 
of ^ and the four parameters. And, we construct five 
ratios of observables which receive a contribution from 
critical fluctuations that is independent of ^ and inde- 
pendent of all four currently poorly known parameters. 
This means that we make robust predictions for these five 
ratios, predictions that could be used to provide a strin- 
gent check on whether enhanced fiuctuations discovered 
in some experimental data set are or are not due to criti- 
cal fluctuations. We close in Section IV with a discussion 



of remaining open questions. 

We shall find that critical fluctuations can easily make 
contributions to the higher moments of the proton multi- 
plicity distribution that are larger than those in a Poisson 
distribution by more than a factor of 100. In Appendix A 
we convince ourselves that we can construct a reasonable 
looking, but somewhat ad hoc, distribution whose higher 
moments are this large. What we are able to calculate in 
Section II is moments of the distribution, not the distri- 
bution itself. In Appendix A we construct a toy model 
distribution that has moments comparable to those we 
calculate. We also use this toy model to obtain a crude 
gauge of how our results would be modified by any effects 
that serve to limit the maximum proton multiplicity in a 
single event. 

In Appendix B we apply our calculation to determine 
the contribution of critical fluctuations to the third and 
fourth cumulants of the event-by-event distribution of the 
mean transverse momentum of the pions in an event. We 
find that the critical contribution to these non-Gaussian 
cumulants are quite small, smaller even than the con- 
tributions of Bosc- Einstein statistics. For this reason, 
throughout the main text of the paper we focus entirely 
on number fluctuations, rather than transverse momen- 
tum fluctuations. 



A. Moments and cumulants of fluctuations 

We expect to see a peak in the Gaussian and non- 
Gaussian cumulants of particle multiplicity distributions 
near the critical point as we change y/s. In this subsec- 
tion, we describe how to calculate these higher cumulants 
from experimental data. 

Consider an ensemble of events in each of which we 
have measured the number of particles of two species, 
which we shall denote x and y. The possibilities for x 
and y that we consider later include the number of pi- 
ons iV^, the number of protons iVp, and the number of 
protons minus antiprotons Np-p = Np — Np. In each 
case, the number that is tallied should be the number 
of particles of the desired species near mid-rapidity in 
a specified window of rapidity. This window in rapid- 
ity should be at least about one unit wide, in order for 
our results to apply without significant acceptance cor- 
rections [25j . Furthermore, the longitudinal expansion of 
the matter produced in the collision reduces correlations 
among particles separated by much more than one unit 
in rapidity [251, making larger windows unnecessary. 

We denote the average value of x and y over the whole 
ensemble of events by {x) and (y) . Throughout this pa- 
per, we use single angle brackets to indicate the ensemble 
average of a quantity whose event-by-event distribution 
has been measured. And, we shall denote the deviation 
of X and y from their mean in a single event by 

Sx = X — (x) 

Sy = y-{y) (i.i) 
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We now define the cumulants of the event-by-event dis- 
tribution of a single observable, say x. The second and 
third cumulants are given by 



^3. ^ {{x'))^{{6xf] 



(1.2) 
(1.3) 



where we have introduced two equivalent notations for 
the cumulants. The second cumulant K2x is the variance 
of the distribution, while the skewness of the distribution 

is given by n-ixji^-ix ■ The fourth cumulant is different 
from the corresponding fourth moment: 



Ai4. = ((a;4)) = ((5:E)4)-3((fe)2)^ 



(1.4) 



The kurtosis of the distribution is given by yi^xl ^^x■ 

The defining property of the cumulants is their addi- 
tivity for independent variables. For example, if a and 
h are two independent random variables, then Ki(a+f)) ~ 
i^ia + Hib- This property is easily seen from the cumulant 
generating function 



5(/i)=log(e''^-) 



(1.5) 



which is manifestly additive. The ri'th cumulant of the 
x-distribution is given by 



9^" 



(1.6) 



/j=0 



Using the double bracket notation introduced above, 
.g(/i) = {{e^^)). As a result of their additivity, cumulants 
of extensive variables, such as Nj, or N^^, are all them- 
selves extensive, meaning that they are proportional to 
the volume of the system V in the thermodynamic limit. 

We shall also consider mixed cumulants, which gener- 
alize the more familiar Gaussian measures of correlations 
to non-Gaussian measures. These are generated by 



dip-, v) 



n.m 



log(e 



^ 5x+y 5y\ 



and, for example, are given by 



Kixiy = {{xy)) = {SxSy) , 
Kix2y = {{xy'^)) = {5x{5yf) , 

I^2x2v = {{X^y^)) 

^{{5xf[5yf)--2{5x5y)^ 



(1.7) 



(1.8) 
(1.9) 



{5xf 



{{xy')) 

{Sx{Sy)'] 



3{SxSy) {{Sy)' 



{5y)') , 

(1.10) 



(1.11) 



For two extensive variables x and y such mixed cumulants 
are also extensive, proportional to V . 

We have described how to obtain the cumulants Kij,, 
Kjy and Kixjy from a data set consisting of an ensemble 
of events in each of which x and y have been measured. 



We can now define the intensive normalized cumulants 
that we shall analyze: 



(1.12) 
(1.13) 
(1.14) 
(1.15) 
(1.16) 



{N.) ' 

{Np + Np) ' 

^ipjn 

{Np + Np)-/r{N^Y/r ' 



where r = i + j. 



If Nt^, Np and Np are statistically independent and 



Gaussian distributed, then the u>2^s in (1.12), (1.13) and 



(1.14) are nonzero and all the other w's vanish. 



If N^, Np and Np are statistically independent and 



Poisson distributed, then all the Wi's in (1.12), (1.13) 



and (1.14 1 with i > 2 are equal to 1, and all the mixed 



cumulants vanish and therefore so do the w's in (1.15) 
and drTel). 



In this paper we shall calculate the contributions of 
critical fluctuations to the normalized cumulants (1.12 1, 



(1.13) and (1.14) for i = 2, 3 and 4 and the normalized 



mixed cumulants (1.15) and (1.16) for i's and j's such 
that r = 2, 3 and 4. 



B. Dependence of ^ on /is 

We shall close this Introduction (in Section I.G) by 
illustrating possible experimental outcomes of measure- 
ments of the cumulants defined in Section LA, assuming 
that the matter produced at the freezeout point of the 
fireball evolution for some collision energy ^/s is near the 
critical point. In Section I.C we shall present only results, 
while the calculations involved are presented in Section 
II. What we shall calculate in Section II is the contribu- 
tion of critical fluctuations to the observables defined in 
Section LA, in terms of the correlation length ^. In order 
to give an example of possible experimental outcomes, we 
need to make an illustrative choice of how the correlation 
length ^ that is achieved in a heavy ion collision depends 
on 

To start, let us assume that the critical point occurs 
at 11% = 400 MeV. Let us also assume that because the 
fireball only spends a finite time in the vicinity of the 
critical point the correlation length reaches a maximum 
value of ^max = 2 fm in the collisions in which the freeze- 
out point is closest to the critical point during an energy 
scan. We stress that our choices of iig and fmax are 
arbitrary, made for illustrative purposes only, and are in 
no way predictions. 

How does the correlation length achieved in a heavy 
ion collision depend on the /i^ at which the matter pro- 
duced in the collision freezes out? Glose to the critical 
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point, the equilibrium correlation length ^eq is very long 
and there is not sufficient time for the actual correlation 
length ^ achieved in a collision to reach ^gq [21] • Lets sup- 
pose that ^ reaches ^cq for — A^bI ^ for some W, 
while for I/is — /i^l ^ ^ finite time effects limit ^ such 
that it peaks at ^max- In principle, ^cqlMs) could one day 
be determined from lattice QCD calculations, but these 
calculations are challenging at 7^ because of the no- 
torious fermion sign problem, so this day remains in the 
future. At present, all we can do is require that the static 
correlation length ^oq satisfy the constraints imposed by 
the universality of critical behavior at long wavelengths. 
The universal behavior is really only attained in the limit 
in which — > and ^max — >■ 00, so our use of it in the 
present context is illustrative but not quantitative. As 
a function of /is ~ Mb' the universal regime ^oq niust 
scale as ^ — /±|/is — where is the relevant 

critical exponent^ and /+ and /_ are the amplitudes of 
the singularity on the crossover and first-order side of 
the transition respectively. The precise value of the criti- 
cal exponent is v = (2 — a) /3 ~ 0.63, with the numerical 
value being that for a critical point in the Ising universal- 
ity class But, in our calculation in Section II we shall 
be neglecting the small anomalous dimensions associated 
with nonvanishing values of the exponents 77 ~ 0.04 and 
a ~ 0.1. So, to be consistent, here too we shall sim- 
ply use 1/ = 2/3. The ratio of the amplitudes f+/f- is 
also a universal quantity. In the Ising universality class, 
/+//_ « 1.9 [27]. Since /+//_ > 1, the correlation 
length falls off more slowly on the crossover side ^ < fig . 

The simplest ansatz for £,{ijlb) that we have found that 
incorporates the physics that we have just described is 



w(tiBY 



1/3 



with 



W{fiB) = W + 6Wta.nh 



fJ-B - fJ-B 



(1.17) 



(1.18) 



^ For our illustrative model of the ^{(j-b) dependence along the 
freezeout curve we are assuming that where the freezeout curve 
passes the critical point it is approximately parallel to the tran- 
sition line (crossover and first-order lines). The region of the 
QCD phase diagram in the (fis , T) plane near the critical point 
can be mapped onto the Ising model phase diagram, whose 
reduced temperature and magnetic field axes are convention- 
ally denoted by t and h, respectively. Upon approaching the 
Ising critical point along the t-direction, i.e., along the tran- 
sition line, ^eq ~ t"" ~ 4"^/'^, while along the /i-direction, 
^eq ~ h-"/P^ ~ As long as h < t''* on the freezeout 

curve, the t-like scaling dominates and, since — ~ t, we 
obtain feq ~ Ia'b ~ Msl""- condition h <C is violated 

at points on the freezeout curve that are very close to the criti- 
cal point, t 0, where the h-like scaling sets in. For simplicity 
we assume that this small-t segment of the freezeout curve in 
the QCD phase diagram lies in a region where the equilibrium 
correlation length ^eq already exceeds ^max = 2 fm, and thus 
^ ~ fmax in this segment. 
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FIG. 1: The correlation length ^(/is) achieved in a heavy 
ion collision that freezes out with a chemical potential /xs, 
according to the ansatz described in the text. We have as- 
sumed that the collisions that freeze out closest to the critical 
point are those that freeze out at fi% = 400 MeV. We have 
assumed that the finite duration of the collision limits ^ to 
i < ^max = 2 fm. We show ^(/is) for three choices of the 
width parameter A, defined in the text. The choices of pa- 
rameters that have gone into this ansatz are arbitrary, made 
for illustrative purposes only. They are not predictions. 



where W and w are nonuniversal parameters to be chosen 
and SW is specified by requiring that 



W 

W 



SW 



3/2 



= 1.93/2 



(1.19) 



We have constructed (1.17) such that ^ has the universal 
behavior of ^eq when \fiB — ^ W{ijLb), but has a 
peak that is cut off at ^ = Cmax where fiB = fJ's- We 
have chosen the shape of ^ in the vicinity of the peak 
arbitrarily, for illustrative purposes, not via analysis of 
the rate of growth of ^ during the finite duration in time 
of a heavy ion collision. In Fig. [T] we show two instances 
of our ansatz for ■^(/is). They differ in their choice of the 
width of the peak. We shall define the width A as the 
distance in /x^ between the two points at which S,{hb) 
crosses 1 fm, i.e. the width in /is within which ^ > 
1 fm. The three curves in the figure have A=50, 100 and 
200 MeV. In all three cases we have chosen w = 0.1 A. 
(With this choice, W = 0.189A and 5W = 0.084A.) 
There is no reason to expect that A should be small and, 
indeed, in model calculations it seems to be larger than 
100 MeV [28 . Ultimately A should be determined by 
lattice calculations; one first attempt to do so indicates 
A ~ 100 MeV [IZllMj- 



C. Cumulants near the critical point 

We shall concentrate our analysis on observables char- 
acterizing the fluctuations of pious and protons. Pions 
are the most abundant species produced in relativistic 
heavy ion collisions. Protons are important, among other 
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FIG. 2: The ^s-dependence of u^p, the normaHzed 4th cu- 
mulant of the proton number distribution defined in (1.131, 
with a /^fl-dependent ^ given by | |1.17[ ). We only include the 
Poisson and critical contributions to the cumulant. In the top 
panel we choose fig = 400 MeV and illustrate how aj4p is af- 
fected if we vary the width A of the peak in ^ from 50 to 100 
to 200 MeV, as in Fig. [l] The inset panel zooms in to show 
how ojip is dominated by the Poisson contribution well below 
fi%. In the lower panel, we take A — 100 MeV and illustrate 
the effects of changing ii% and of reducing the sigma-proton 
coupling (?p from our benchmark gp = 7 to Qp — 5. 



reasons, because their fluctuations are proxy to the fluc- 
tuations of the conserved baryon number |30j and be- 
cause their coupling to the critical mode a is relatively 
large. 

We have defined the normalized cumulants of the pro- 
ton and pion distributions in ( |1.13 1 and (1.12 1 and the 
normalized mixed cumulants in (1.15). Fig^ shows how 



LOip might look like, with ^(/^s) given by Eq. (1.17 1. We 
illustrate how uj-^p changes if we vary the location of the 
critical point and the width A of the peak in Fig. [TJ 
as well as the sigma-proton coupling gp. As we shall see 
in Section II A there are four nonuniversal parameters 
that (for a given ^max) govern the height of the peaks 
of the normalized cumulants. These include gp and the 
sigma-pion coupling G, as well as two p arameters A3 and 
A4 that we shall define in Section II A We have used as 
our benchmark values G ~ 300 MeV, g — 7, X3 = 4 and 
A4 = 12. As we shall discover in Section II and discuss 
at length in Section III, the heights of the peaks of dif- 
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FIG. 3: The ^_g-dependence of selected normalized cu- 
mulants, defined in 11.121, (1.131 and Jl. 15 1, with a /is- 
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dependent ^ given by (1.171 as in Fig. Ill We only include 
the Poisson and critical contributions to the cumulants. We 
have set all parameters to their benchmark values, described 
in the text, and we have chosen the width of the peak in 
^ to be A = 100 MeV. Note the different vertical scales in 
these figures and in Fig. [2] The magnitude of the effect of 
critical fluctuations on different normalized cumulants differs 
considerably, as we shall discuss in Sections II and 111. As we 
shall also discuss in those Sections, ratios of the magnitudes 
of these different observables depend on (and can be used to 
constrain) the correlation length ^, the proton number den- 
sity Hp, and four non-universal parameters. We shall also see 
in Section III that there are ratios among these observables 
that are independent of all of these variables, meaning that 
we can predict them reliably. For example, we shall see that 
critical fluctuations must yield 1^2^271 ~ {^4p — l)(t<J47r — 1) and 

t^iplTT = i^3p - l)^(t^37r - 1) and Ujfp2„ = {i^3p - l)('^37r " 1)^. 

(The subtractions of 1 are intended to remove the Poisson 
background; in an analysis of experimental data these sub- 
tractions could be done by subtracting the ujip or ujj-n de- 
termined from a sample of mixed events, as this would also 
subtract various other small background effects.) 



ferent normalized cumulants are affected differently by 
variations in these four parameters. Fig. [3] shows how six 
more different normalized cumulants vary with hb- In 
this figure we keep all parameters set at their benchmark 
values, deferring a discussion of how these peaks change 
with parameters to Section III. 

In the case of free particles in the classical Boltzmann 
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regime, with no critical fluctuations, the fluctuations of 
any particle number obey Poisson statistics. The Pois- 
son contribution to w^p and w^jr is 1, and in the figures we 
have added this Poisson contribution to the contribution 
from critical flucuations that wc calculate in Section II. 
There is no Poisson contribution to the mixed cumulants 
ijJipjTi- In reality, in the absence of any critical fluctua- 
tions the 1 of Poisson statistics gets few percent contri- 
butions from Bose-Einstein statistics, from initial state 
correlations that are incompletely washed out, and from 
interactions other than those with the fluctuations that 
are enhanced near the critical point. System size fluctua- 
tions are also a potential non-critical contribution to the 
fluctuation measures. We do not attempt to estimate 
this effect (see, e.g., Refs. [5D], [21]), and assume that 
sufficiently tight centrality binning suppresses it. We are 
ignoring all of these non-critical corrections to the Pois- 
sonian 1 and in the plots shown here we only include the 
Poisson and critical contributions to the cumulants. Ex- 
isting data on K,i[p-p) / ii2{p-p) — 19.6, 62.4 and 
200 GeV |32] confirm that the non-critical corrections to 
the Poissonian 1 are indeed small, and confirm that it is 
possible to measure 4th order cumulants with an error 
bar that is much smaller than 1. 

We can clearly see the peak in all the normalized cu- 
mulants near the critical point. In many cases, the peak 
due to critical fiuctuations is larger than the Poisson con- 
tribution by more than an order of magnitude.^ The 
results indicate that the more protons are involved in 
the observation measure, the easier it is to identify the 
critical contribution. The reader who would like to see 
an example of a probability distribution that has UJ4 as 
large as uJ4p gets in Fig. [2] should consult Appendix A. 
A more comprehensive discussion of the results is given 
in Sections |lll] and |IV[ but it is readily apparent that 
the measurement of these observables in heavy ion col- 
lisions at a series of collision energies is very well suited 
to ruling out (or discovering) the presence of the QCD 
critical point in the vicinity of the freeze-out points of 
the collisions in such an energy scan. 



II. CALCULATING CRITICAL CORRELATORS 
AND CUMULANTS 

In this section, we show how to calculate the critical 
point contribution to the cumulants of the particle mul- 
tiplicity distribution of pious, protons and net protons. 
We essentially show how to obtain the normalized cumu- 
lants in Figs. [2] and [3] as the location of the critical point, 
fig, changes. We begin in Section II A by calculating the 



^ Although it is a small effect, note that the peaks of any of the 
cumulants involving protons do not occur exactly at the fig at 
which ^{ij,b) from Fig. [l] peaks, because the cumulants them- 
selves depend directly on the proton number density and hence 
on fig, as we shall see in Section II. 



correlators that describe the critical contributions to the 
fiuctuations of the occupation number of pions and pro- 
tons with specified momenta. We use these correlators 
to calculate the normalized cumulants in Section III Bl 



A. Critical point contribution to correlators 

Fluctuations of observables, such as particle multiplic- 
ities, are sensitive to the proximity of the critical point if 
the particles under consideration interact with the crit- 
ical field a — the field whose equilibrium correlation 
length diverges at the critical point. In this Section, 
we shall treat the a correlation length ^ as a parame- 
ter, in this way avoiding any consequences of our lack of 
knowledge of the dynamics of how the long wavelength 
correlations in the a field grow. In order to use the re- 
sults of this section to make the plots in Section I.C, in 
Section LB we had to make an ansatz for S,{hb)- But, 
the results of this section, expressed in terms of ^, are 
independent of the uncertainties in that ansatz. 

We can describe the fiuctuations of the cr-field by a 
probability distribution of the form 



P(a)^exp(-r!((7)/T), 



(2.1) 



where is the effective action functional for a. It can be 
expanded in gradients and powers of a as 



{Vaf 



A 



A4 



(2.2) 



2' ' 2 3 4 

In this expression the sigma-field screening mass is 

= (2.3) 



and, near the critical point, the and interaction 
couplings are given by 

A3 = A3r(T0-3/2^ and A4 = A4(T0"\ (2.4) 

where the dimensionless couplings A3 and A4 do not de- 
pend on ^, but do depend on the direction of approach 
to the critical point, as described in Ref. [21]. These cou- 
plings (and their dependence on direction) are universal 
and they have been determined for the Ising universality 
class [33]. Throughout this paper we shall use A3 = 4 
and A4 = 12 as benchmark values, because these are the 
midpoints of the ranges of values known for these con- 
stants [211 [33] ■ In fact, both A3 and A4 will vary with 
/is, as the location of the freeze-out point moves in the 
phase diagram, relative to the critical point. We shall 
not attempt to parametrize the /is-dependence of these 
parameters, however, because the dominant source of /in- 
dependence in our results is the variation of ^ with /is, 
and our knowledge of ^(/is) is sufficiently uncertain (as 
we saw in Section LB) that this uncertainty would domi- 
nate any increase in precision that would be obtained by 
modelling the /is-dependence of A3 and A4. 
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The correlation functions and fluctuation moments and 
cumulants of the critical field a itself can be calculated 



directly using the probability distribution given in (2.11, 
but these quantities are not directly observable. The 
long wavelength fluctuations in the a-field manifest them- 
selves in observable quantities in so far as they affect the 
fluctuations of the occupation numbers of particles that 
couple to the cr-field. This coupling to the fluctuating 
field a contributes to the moments of particle fluctuations 
the terms proportional to the corresponding moments of 
a itself [23]. Both protons and pions couple to the a 
field. We shall define the strengths of the corresponding 
couplings Qp and G through the respective terms of the 
effective Lagrangian (following the notations of [20 l l30 j ): 



gp a pp. 



(2.5) 



where tt^ is the (charged) pion field and p is the Dirac 
fermion field of the protons. The coupling that we denote 
gp is often just called g. We shall make the discussion 
that follows similar for protons and pions by defining a 
dimensionless measure of the sigma-pion coupling 



G/tTIt, 



and using the notation g when we intend an equation to 
be valid for either pions, with g — !■ g^^, or protons, with 
g ^ gp- Throughout this paper we will use G — 300 
MeV (see Rcf. [20_ for a discussion of how to estimate 
G) and gp — 7 (see, e.g., j34|) as benchmark values. It 
is important to bear in mind that both these parame- 
ters and A3 and A4 are all uncertain at the factor of 2 
level. These parameters enter into our calculations of 
the various normalized multiplicity cumulants, making 
absolute predictions of these observables in terms of ^ 
difficult. The advantage that we have, however, is that 
we will be able to calculate many different normalized 
cumulants that depend differently on these parameters. 
In Section |III| we shall discuss how to use deliberately 
chosen ratios of cumulants to measure and even overcon- 
strain various combinations of these parameters. And, 
we shall find five ratios of cumulants that are indepen- 
dent of the values of all of these parameters, allowing us 
to make parameter-free predictions of these ratios. 

The critical contribution to the proton or pion corre- 
lators arises from virtual a-exchanges which introduce 



powers of the correlation length ^ 



where 



is the cr-field screening mass. As the correlation length 
grows in the vicinity of the critical point, the contribution 
to the particle correlators due to a cr-exchange dominates 
over other non-critical contributions. The effect of such 
an interaction on the two-point particle correlators was 
studied in Refs. [201 121] and on higher-point correlators 
in Ref. [H]. In this subsection we will only look at the 
particle correlators and in the subsequent sections we will 
show how to calculate cumulants of particle multiplicity 
distributions from the correlators. 

The contribution of critical fluctuations to the 2-, 3- 
and 4- particle correlators due to cr-exchanges can be 



calculated using the diagrammatic method developed in 
Ref. ^ (see also Refs. ^20] and ^). We shaU write 
the correlators using a notation that applies to either 
protons or pions. They describe the correlation between 
the dn^s at different momenta, where Sn^ = n^ — (rik) is 
the difference between the occupation number of the k'th 
pion or proton mode in momentum space in a particular 
event and its mean value. The correlators are given by 



9' 



ki 



mlV T 

VT 7k2 



(2.7) 



2(PX 
2d^X3 

1/2^3/2 



3 / 9_\^ "ki ^k2 ^k 
2 



m^a J 7ki 7k2 7k3 

^3^9/2< <<^(2., 



(2.6) {{Sn]^^Snk^Sni^Jnk^)), 





7ki 7k2 7k3 



4 9 9 9'? 

7ki 7k2 7k3 7k4 



. 2A^ - A4 1 gr ^ , (2.9) 

y^T^ ^ ' 7ki 7k2 7k3 7k4 



where we have used (2.3) and (2.4 1 and where we must 



now explain many aspects of our notation. The subscript 
CT indicates that we have only calculated the contribu- 
tion of the critical fluctuations to the correlators. The 
double brackets around the quartic correlator indicate 



that what is evaluated is the cumulant, as in (1.4). The 



equations (2.7), (2.8), (2.9) apply to both protons (with 
g — gp) and pions (with g — g-n — G/to^). The degener- 
acy factor d is 2 for both protons and pions. (For protons, 
d — dp — 2 counts the number of spin states. For pions, 
d — d^^ — 2 counts the number of charge states — tt"*" 
and TT^. These degeneracy factors appear because the 
coupling to the cr-field is both spin and charge "blind".) 
The variance of the fluctuating occupation number dis- 
tribution is denoted by and is given by 



i;2 = (nk)(l±(nk)). 



where, as usual, 



exp[(7km-^)/T] =F 1 



(2.10) 



(2.11) 



with m = m^, /i = and the upper sign for pions and 
m = nip, fi = and the lower sign for protons. And, 
finally. 



7k 



\/k2 



(2.12) 
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is the relativistic gamma-factor of the particle with mas s 
m with a given momentum k.'^ We see from Eqs. (2.7)- 



( 2.9 ) that these correlators, and hence the cumulants that 



we will obtain from them, are proportional to powers of 
the correlation length ^ and so peak at the critical point. 

Now let us turn to mixed pion-proton correlators. The 
2 pion - 2 proton correlator is given by 



Now let us focus on how one can obtain higher cumu- 
lants of the particle multiplicity distributions using the 
correlators found in the previous section. As an example, 
let us evaluate the critical contribution to the normalized 
fourth cumulant of the proton multiplicity distribution, 
W4p defined in (1.13). The total multiplicity Np is just 
the sum of all occupation numbers rik, thus (see ref. |24j ) 



{{iSNpr))^ 



(2.14) 




2 TT 2 77 2 P 2 p 2 

^ki ^k2 ^k3 "k4 

7ki 7^, 7k3 7k4 



= V' 



k2 •'ks Jk4 



(^12 l\ 2 2.7 <1 <^ <l «k4 (r. . ON 

^ (2A3 - j (2-13) 

The prescription for obtaining other mixed correlators 
from the correlators (2.7 - 2.9) should be clear: each As we discussed in Section LA, see (1.131, we normalize 
particle brings its own corr espo nding factor dgv^/j]i to 
the expression in, e.g., Eq. (2.9). In this way, the 1 pion 
- 3 proton and 3 pion 
be obtained from Eq. 
pion 



where 



d^k 
(27r)3' 



(2.15) 



(2.9) 



1 proton mixed correlators can 
the 1 pion 



2 proton and 2 



1 proton mixed correlators can be obtained from 
1 proton can be obtained from 



the cumulant by dividing by the total proton multiplic- 
ity Np. To simplify notation below, it is convenient to 
introduce the proton and pion number densities 



(Np) 



Eq. (2.8), and the 1 pion 
Eq. (2.7). 



Another useful fluctuating quantity to consider is the 
net proton number correlators (the net proton number is 
defined as the number of protons minus the number of 
anti-protons: Np-p = Np — Np). In order to obtain the 
corresponding correlators one can begin with the similar 
correlators for the protons and replace u^^ with (w^^ — 
^k^)' where is the occupation number variance for 
anti-protons. (See, e.g., Ref. [50]). 

In the next section we will use these correlators to eval- 
uate cumulants of particle multiplicity distributions for 
pions, protons and net protons and see how they can be 
used to locate the critical point. 



V 



V 



I \"-k/ 

/k 



TT 



^ /"OO 

~2 



dEEjE^ 



d-^ 

k 

1 f'^^dEE^/W 



TT 



,E/T _ I 



(2.16) 



(2.17) 



The result we find for the normalized cumulant can then 
be written as 



6(2Ai-A4) ,, ( fvf; 



(2.18) 



B. Energy dependence of pion, proton, net proton, 
and mixed pion/proton multiplicity cumulants 

In this section we will concentrate on cumulants of the 
particle multiplicity distributions and how they vary as 
we change the location of the critical point and change 
the value of parameters. Another application of the cor- 
relators given in the previous section is the calculation 
of the critical point effect on higher moments of the fiuc- 
tuation of mean transverse momentum pT- We find that 
the critical contribution to pT fluctuations is rather small 
(e.g., smaller than the enhancement due to Bose statis- 
tics) and thus not as useful in the search of the critical 
point. Details can be found in Appendix [B| 



^ A note on subscript/superscript notation: we denote momentum 
subscripts with a bold letter k. Subscripts/superscripts denoting 
particle type, e.g. p for protons, will be in normal typeface. 



We can see from expressions (2.7) - (2.9) that higher cu- 



mulants are proportional to higher powers of and thus 
increase by a larger factor near the critical point where ^ 
becomes large. For example, the third and fourth cu- 
mulants are proportional to ^^/■^ and respectively. 
If the correlation length ^ increases from ~ 0.5 fm to 
Cmax = 2 fm as in Section LB, these cumulants are sub- 
stantially enhanced — as we have seen in the plots in 
Section I.C. 

With an explicit expression for w^p^a in hand, we can 
now write our general result for uJipjTr,a 
can also include uj., 

( 1.12 ) in the notation via setting j = or i 



in (1.15). We 



jp^cr and WjTr.o- deflned as in (1.13) and 



in oji. 



We obtain 



to. 



ipjTT 



'Ji,0 ■ 



^3,0 



A;,(r-1)! a. 



prcfactor 
ipj-TT 



— £2 
i/r j/r ^ 



(2.19) 



no 




FIG. 4: Proton number density Up and net proton number 
density at chemical freezeout as functions of 

fiB- Both depend on T as well as fiB', we have taken T{iib) 
as in (|2?24l). 



We have normalized both rip and np_p u sing the 



constant no of (2.231 introduced in (2.191 and (2.201 



where we have defined 



prcf; 



'actor 1 



}r-3 



i/r j/r 
Tin TZtt 



and 



(2.20) 



,7r 2 



/k 7£ 



dp9p I (2.21) 
/k 7k 



A2 = 1, A3 = A3 and a; = 2A^ - A4. (2.22) 



In the second line of (2.19) we have factored out the two 
main sources of fiB dependence: the correlation length ^ 
depends on /i^ as we have discussed at length in Section 
IB and, if the normalized cumulant involves the proton 
multiplicity it depends on Hp, which increases rapidly 
with increasing /i^ as shown in Fig. |4] We have denoted 
all of the remaining factors in our result for the contribu- 
tion of critical fluctuations to the normalized cumulant 
by 1^'"^*°'^ , which depends only weakly on hb as we 
illustrate in Fig. [Sj The number density no is an ar- 
bitrary constant — note that it cancels when (2.20) is 
substituted into ( |2.19[ ) — introduced in order to make 
i^f^fj"^*'"'^ dimensionless. We shall choose 



no 



(5 fm)"^ 



6.116 X lO-'^ GeV^ 



(2.23) 



With this choice, {np)/no is of order 1 at the /is of in- 
terest to us — see Fig. |4] — and none of the different 
u}fpj^^'^*'°'^s are orders of magnitude smaller or larger than 
1, as illustrated in Fig. [5] 

Let us now walk through the physics behind the dif- 
ferent pieces of the expression (2.19). The Kronecker 



deltas describe Poisson fluctuations, which are of course 
^-independent. As we described in Section lA, they con- 
tribute 1 to the Wip's and the watt's and they make no 



prefactor 

1.6 f I 

1.4 " 2;r 

IplTT, l(p-p)l;r 

0-8 : __ 2p,2(p-p) 

0.6 I 

0.4^ _____ 

' /ifiVGeV 

0.0 0.1 0.2 0.3 0.4 0.5 ^ ' 

prefactor 
^^3 



371 

lp2x, l(.p—p)2n 
3p, 3(p-p) 




4n 

lp3;r, Hp-p)3n 
2p2x, 2(p-'p)2jr 
3pl;r, 3(p-'p)ln 
4p, 4(p-p) 



FIG. 5: The Ms-dependence of Lof;f:^'°' and de- 
fined in ( |2.19[ ), \2.20\ and ( |2.25[ ). The three panels are for the 
normalized cumulants with r = i+ j = 2, 3 and 4, respec- 
tively. The curves can be used to determine how the height 
of the peak in the critical contribution to the normalized cu- 
mulants changes as we vary /x^, the /xs at which ^ = fmax 
and at which (to a very good approximation) the normalized 
cumulant has its peak. The height of the peak in cj^pj^ [or 



is proportional to (np/no)'"*' [or {up-p/no) 
multiplied by the prefactor plotted in this Figure. We have 
taken T{iib) as in (2.241 and have used the benchmark pa- 
rameters G = 300 MeV, gp = 7, A3 = 4 and A4 = 12. 



contribution to the mixed cumulants in which i and j are 
both nonzero. More realistically, the 1 of Poisson statis- 
tics gets few percent contributions from Bose-Einstein 
statistics (which are calculable), from initial state corre- 
lations that are incompletely washed out, and from inter- 
actions other than those with the critical u-mode. We are 
ignoring all of these noncritical corrections to the 1. In 
principle, with sufficiently precise data their magnitude 
could be measured far away from the critical point and 
this background could then be subtracted. If this back- 
ground were significant, one could also try to study and 
calculate these corrections theoretically. Present data on 
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Ki{p-p) / K2{p-p) at y/s = 19.6, 62.4 and 200 GeV indicate 
that the corrections to the Poissonian 1 are quite small, 
but this should be investigated also for othe r cumulants. 

The second, ^-dependent, term in (2.191 is the con- 
tribution to uJipj.,^ made by the critical fluctuations. It 
grows proportional to near the critical point. 

We see evidence of this in the heights of the peaks in dif- 
ferent w's in Fig. [3j but it is also clear from this Figure 
that the r-dependent difference in the power of ^ is not 
the on ly ini portant source of /i-dependence. Indeed, we 



are obtained is the replacement of ^ with vZ^ 



„P2 



We find 



r I r . prefactor / '^p—p 



no 

|r-3 



(2.25) 



see in (2.20) that i^f^f^'^*' is proportional to Up and to 



ttp and, it turns out, ap/up is close to constant. This 
means that the dominant /is-dependence of the critical 
contribution to iOipj^^ at a given ^ is Up which we 



have factored out in (2.19) making the /^s-dependence 



in ujf^f^^^ rather mild. We can see the 



-i/r 



where Up-p = Up — np is the net proton number den- 
sity. In comparison with Eq. (2.19), we have pulled out 

another factor, ( "V'' I , which describes the vanish- 

ing of the critical contribution to net proton cumulants 
at = 0, see Fig. |4j It then turn s out that the pref- 
actor uj^'^^_^-y'^ (defined as in t.20\, but muhiplied by 



"p+Kp \ ^j^j with wP^ replaced by v 
from w; 



p2 
k 



vV) differs 



prefactor 



jpj^ by less than half of one percent, which is 
less than the thickness of the curves in Fig. [5j Hence, 

, 1 1 J ■ 1. prefactor 

these curves also depict lo , . . 

f l{p-p)]TT 



■p — de- 
pendence of the height of the peak in uj^p in the lower 
panel of Fig. [2] in this figure ^max is the same for all the 
curves so the /^s-dependence of the height of the peaks 
in uj4p comes from its rip-dependence. 

For i = 0, meaning for a cumulant involving pions 
only, there is no large np-dependence in cj^^ and the In order to evaluate either ( |2.19[ ) [or ( |2.25[ )] and com- 
height of the peak in a figure like Fig.[3]is proportional to 
^prefactor ^ and the dominant /is-dependence oiujjT^ itself 
comes from its dependence. For observables involv- 
ing protons (z > 0) , the dominant contribution to the /i^ 
dependence of the height of the peak in u comes from the 
factor Up , and the slowly varying prefactor in Fig. [s] 
adds relatively little to that strong dependence. 

In plotting the curves in Fig. |4] and Fig. [5j we have 
allowed for the fact that the chemical freeze-out tem- 
perature T decreases somewhat with increasing fj,B- 
We have described this dependence using an empirical 
parametrization of heavy ion collision data from Ref. |35j : 



pare to data, we need the proton number density Up [net 
proton number density Up-p] at each collision energy ^/s. 
These can be extracted from data via the conventional 
statistical model fits done at each -^s that give /is and 
T at chemical freeze-out at each ^/s. The value of Up at 
chemi cal freeze-out is specified in terms of and T by 
( 2.16 ) and the value of np is given by the same expression 
with fiB replaced by —fiB, so these number densities can 
also be obtained from data. So, at each collision energy, 
one should take the /is and T from the statistical model 
fit, evaluate rip and np_p, and then plug these into (2.19) 
and (2.25) and see what conclusions can be drawn about 



^(ms) 



a 



(2.24) 



with a = 0.166 GeV, b = 0.139 GcV"! and c = 
0.053 GeV~'^. Almost all of the /is-dependence of the 
i^fpf^^°'^s plotted in Fig. actually comes from the fiB- 
dependence of the chemical freeze-out temperature T. In 
plotting Fig.[5j we have used our benchmark values of the 
four nonuniversal parameters that determine the uJipj^^ 
for a given ^, namely g^r — G/rriTr = 2.1, gp = 7, A3 = 4 
and A4 — 12. 

Finally, completing our discussion of the proton-pion 
cumulants and Fig. [5] we note that in the lower panel in 
this figure there is a point where all five ujfpj^^^°'^s with 
r = 4 cross. This occurs if at some value of hb it so 
happens that aTr{no/n^y^'^ and apUo/rip coincide. 

We now turn to the net proton multiplicity distribu- 
tion, where by net protons we mean Np_p = Np — Np. 
The calculation of the normalized cumulants involving 
the net proton multiplicity and the pion multiplicity, 
namely (1.16), is analogous to the calculation we have 
described above. As we discussed in Section IIA, the 
only change in the correlator from which the cumulants 



^ and the constants gp, 5^, A3 and A4 using data on as 
many of the normalized cumulants ujipjj^ and uji(^p_p)jT^ 
as possible. We shall provide tuned strategies for this 
analysis in Section |III[ We close this Section with two 
straightforward observations. 

First, the proton/pion normalized cumulant ujipj^^ is 
more sensitive to critical fluctuations than the net- 
proton/pion normalized cumulant i-Oi(^p_p^j„, for any i ^ 
and for any j. As an example let us consider cu^p and 
W4(p_p). We can estimate ij-b{^/s) using the parametriza- 
tion of statistical model fits to data in Ref. [551: 



AiB(\A) = T— 7=, 

1 -I- eys 



(2.26) 



with d = 1.308 GeV and e = 0.273 GeV'^. The proton 
number density np{fiB) is then shown in Fig. H] Then, 
at any fixed value of the correlation length 4 the 7i„- 



dependence that enters the expressions (2.19) and (2.25) 
for cj4p and W4(p_p) 



is 



^ 1 = 0.34, 0.77, 4.9, 31 

no/ 



(2.27) 
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and 



p—p 
no 



i^p—p 



= 0.00072, 0.064, 3.4, 30, 



respectively, when evaluated ai ^ s — 200, 62. 



(2.28) 
19 and 7.7 



GeV. Since Up-p is less than Up (and consequently less 
than Up + Up) at all ^ib — see Fig. |4] — the critical con- 
tribution to LJip is greater than the critical contribution 
to a;4(p_p) at all /is. The analogous argument applies in 



comparing any LUi(^p_p)j^ 



to the corresponding to. 



ipjTT- 



In 



all cases the suppression of the critical contribution to 
Lo^p_p) is most accute at small fj,B, meaning at large y^s. 
We shall ignore a;i(p_p)j^ 



III 



in Section 

Second, we can ask which observable is most sensitive 
to critical fluctuations. For a given ^, the critical contri- 
bution to (jJinjTT is largest when r is largest, since as we see 



from (2.19) this gives w the strongest ^-dependence. The 

demon- 



benchmark values, a^/up = a^/n^ at /i^ — 135 MeV.) 
Both ujip and a;4^ should be measured, and we shall see 
in Section III that if critical fluctuations are discovered 
it will be very important to have data on as many of the 
(jJipjTt as possible. 



III. RATIOS OF CUMULANTS 

In the previous Section, we presented numerical re- 
sults for the contribution made by critical fluctuations to 
various cumulants of particle multiplicity distributions. 
In order to locate the critical point, experimental re- 
sults on multiplicity cumulants will need to be compared 
to the theoretical predictions of the critical contribution 
to these cumulants. But, recall that we had to choose 



where 



experimental measurements reported in Ref. [32] 
strate that K^p_^p-^ / K2{p~p) can be measured with error 
bars that are much smaller than 1, and we expect that 
ojip and 014^ can be measured with comparably small 
error bars. The error bars on measurements of cumu- 
lants with r > 4 will be larger, so until experimentalists 
demonstrate that they can be measured we have focussed 
on cumulants with r < 4. For a given ^ and r, the critical 
contribution to ujipj-^ is largest for z = r if oTpjup > /n^ cumulant, 
or for i = if d^pjnp < a^/riTr- 

It is apparent from Figs. [2] and l3| that L04p ^ Ui-^ at 
= 400 MeV with our benchmark values of gp and g-^, 
meaning that a;4p is the normalized cumulant with r — 4 
that is most sensitive to critical fluctuations. And, it is 
sensitive indeed: we see from the plots in Fig. [2] that if 
^ reaches 2 fm, the critical contribution to co^p will be 
dramatic. Correspondingly, if for example experimental 
measurements were to show that ui^p — 1 < 1 at some 
fiB around 400 MeV, then ^ < 1 fm at that ^b- How- 
ever, if fig is much less than 400 MeV and/or if gp/g-^ 
is much smaller than with our benchmark values, then 
a^p/up could become less than o^^ItIt,, making a;4^ the 
best observable with which to find evidence for the pres- 
ence of critical fluctuations. (With gp and set to their 



benchmark values for four parameters: gp, gT^, A3 and 
A4. These parameters are not known reliably or accu- 
rately enough to permit a quantitative prediction for the 
effect of the critical point on any one of the cumulants we 
have described. The normalized cumulants depend on ^, 
of course, but their dependence on the four poorly known 
parameters would make it difficult to determine ^ from 
data on any one of the cumulants, in isolation. In this 
Section, we suppose that at some y/s there is experimen- 
tal data showing several of the cumulants significantly 
exceeding their Poisson values. We ask how ratios of cu- 
mulants can be used to extract information on ^ and the 
values of the four parameters. And, we construct ratios 
of cumulants that are independent of ^ and all the pa- 
rameters, allowing us to make robust predictions for the 
contribution of critical fluctuations to these ratios. 

The contributions of critical fluctuations to different 
correlators depend on different combinations of ^ and the 
four parameters. For example. 



f^2p,a 
^3p,iT 
f^4p,a 



Vnlgie, 
Vnlgl\^e'\ 

vn^gyx'.e, 



(3.1) 



2~Xl 



A4. For the most general pion-proton 



Vn;g;gi~XU'^'- 



(3.2) 



We 



with r = i + j and with A^ as defined in (2.22). 
have kept the rip-dependence since it introduces signifi- 
cant /is-dependence, but we have suppressed the T- and 
riTT-dependence. In Table |l] we present the parameter 
dependence of various cumulant ratios. Except for the 
first 3 entries, AT^., Np and KipjT^, the quantities we con- 
sider are all ^-independent (i.e. intensive) by construc- 
tion. (In constructing intensive ratios, we can always 
remove V^-dependence by dividing by A'^ to the appro- 
priate power.) Note that although we have not written 
the cr subscripts in the table, the table only describes the 
parameter-dependence of the contributions from critical 
fluctuations. When the ratios in the table are constructed 
from data, the Poisson contribution must be subtracted 
from each measured k separately, before taking a ratio. 
This means that this table will only be useful in the anal- 
ysis of data at values of ^/s at which several k's are differ- 
ent from their Poisson values by amounts large compared 
to the experimental statistical and systematic error bars. 

Looking at Table |lj one can see how to use cumulant 
ratios in order to constrain ^ and the four parameters. 
The correlation length ^ and the four nonuniversal pa- 
rameters always appear in certain combinations in the 
multiplicity cumulants and it turns out that we can only 
constrain four independent combinations. We have con- 
structed the table to highlight ratios that can be used to 
constrain one example of four such combinations, with 
each block delineated by double horizontal lines corre- 
sponding to constraining 
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TABLE I: Parameter dependence of the contribution of crit- 
ical fluctuations to various particle multiplicity cumulant ra- 
tios. We have subtracted the Poisson contribution from each 
cumulant before taking the ratio. The table shows the power 
at which the parameters enter in each case. We only consid- 
ered cases with r = i+j = 2, 3, 4. We defined 2 A3 — A4 = A4. 



ratio 


V Up^flB) gp Qtv A3 Ai 5 




1 _ _____ 


Np 


11 _____ 




I i i j Sr,3 (5r,4 f r — 3 




- I - ^ i j 5r,3 5r,4 \r - 3 


l^ipjTT -^77 l^p 


i j 5r,3 5r,4 \r - 3 


I^2p2itNt7 / KA-nK2p 


- -2 - - -2 


K4pNj^ / K4-jrK2p 


_ -4 _ _ -4 


K2p2TvNp/K4pN^ 


-2 2 - - 


K3pl-^Np/K4pNTr 


-1 1 - - 


,^3/2 , 9/4 ,,1/4 
K3pNp' /K^'p N^' 


-3/2 - 1 - 


l^2pf^4p / l^3p 


- -2 1 


3/2 , 3/2 
'^3p'^27r /^SirKjp 




2 / 2 
^4p/^27r/ ^4nf^2p 




l^4pl^3Tr/ I^4itI^3p 




2 1 
^2p2'K I ^4n f^4p 




3 / 2 
K2plir/ f^3pl^3TT 





1- 57rC — using, e.g., K2p2Tr / KAtt K2p Or 

2- g-^/gp — using, e.g., K2p27vN^ / K4pN^ or 

Np/K4pN^. 

3- ^/9p — using, e.g., K^pNp^^ / kI^^^ N^^'^ . 

4. XyXj — using, e.g., K2pK4p/Kip- 

Since four independent combinations of ^ and the four 
parameters can be constrained by data on these ratios, 
we could, for example, use data to express ^ and three of 
the parameters in terms of the fourth, say gp. 

We can see from Table U that there are also some com- 
binations (e.g., the last five entries in the table) that are 
parameter-independent. The first two of these are in fact 
the ratios of the skewness and kurtosis of protons to pi- 
ous, where skewness and kurtosis are defined as usual 
as 

skewness = ^ , kurtosis = — (3-3) 



* The ratio u)2n = >^2-n I^t: could also be used here. However, we 
have seen that the critical contribution to this quantity is small 
and, given the multitude of alternative choices, we can afford not 
to use this quadratic moment. 



The next row in the table is the ratio of the two rows 
above it, giving a combination that has the virtue that 
it only involves 3rd and 4th cumulants, which is advan- 
tageous since the contribution of the critical fluctuations 
is larger at larger r. The last two ratios in the table are 
quite different, as they involve mixed cumulants, but they 
too are parameter-independent. So, the last five ratios 
in the table have no ^-dependence, no dependence on the 
four poorly known parameters, and no rip-dependence. 
This means that, after we subtract the Poisson contri- 
bution to each of the cumulants involved, we can make 
a robust prediction for the ratios of the contributions of 
critical fiuctuations. We find that these five ratios are all 
precisely 1. 

Now let us see how we can use these five ratios in order 
to locate the critical point. Suppose that as you change 
the center of mass energy -^s of the collisions there is a 
point where many cumulants exceed their Poisson values 
by statistically significant amounts. As we see from Figs. 
[2] and [Sj the qualitative signature of the critical point is 
peaks in the multiplicity cumulants as a function of -^/s. 
Suppose experimental evidence for such peaks begins to 
emerge. The specific ratios of the heights of the peaks 
in Figs. [2] and [3] depended on the benchmark choices for 
parameters that we made in those figures. So, how do 
you check in a parameter-independent fashion whether 
the behavior seen in experimental data is consistent with 
the hypothesis that it is due to critical fluctuations? You 
first subtract the Poisson contributions,^ and then con- 
struct the last flve ratios in Table 1. If the fluctuations 
seen in this hypothetical data are in fact due to the prox- 
imity of the critical point, all five of these ratios will be 
equal to 1, with no theoretical uncertainties arising from 
uncertainty in the values of the parameters. This would 
be strong evidence indeed for the discovery of the QCD 
critical point. 



IV. DISCUSSION 

We have explored the effects of the long wavelength 
fiuctuations that arise in heavy ion collisions that freeze- 
out near the QCD critical point on higher cumulants of 
particle multiplicities. The characteristic signature is the 
non-monotonic behavior of such observables as a function 
of the collision energy, as the freezeout point approaches 
and then passes the critical point in the QCD phase di- 



This can be done by subtracting the values of the cumulants in 
a sample of "mixed events," i.e. events constructed artificially 
from tracks drawn at random from many different events in order 
to remove all correlations. In this way, in addition to subtracting 
Poisson fiuctuations one will also subtract spurious experimental 
effects. The mixed event technique is widely used in the study of 
quadratic moments and it could be used here too, even though 
present data indicate that spurious experimental effects are quite 
small [32]. 
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agram. In Section I we illustrated one possibility for 
how seven different cumulants (third and fourth cumu- 
lants of protons and pions plus three mixed cumulants) 
might behave as a function of /is , the chemical potential 
at freezeout which is the quantity that a beam energy 
scan scans. In Section II, after calculating 21 different 
cumulants as a function of parameters and as a function 
of the correlation length ^ at freezeout, we determined 
that either oj^p or uj^t^ is the most sensitive to critical 
fluctuations, depending on the values of parameters and 
depending on the location of the critical point. However, 
if critical fluctuations are discovered it will be important 
to have measured as many of the cumulant observables 
as possible. In Section III we constructed ratios of ob- 
servables that will allow an overconstrained experimen- 
tal determination of currently poorly known parameters. 
And, we constructed other ratios of observables that, if 
the measured cumulants are indeed dominated by crit- 
ical fluctuations, are independent of all the parameters 
in our calculation and are independent of the value of 
the correlation length. We are therefore able to make 
parameter-independent predictions for these ratios, pre- 
dictions that we hope will some day make it possible to 
determine with confidence that observed fluctuations do 
indeed indicate proximity to the critical point. 

There are several effects that require further investiga- 
tion: 

• In our treatment of the pions we have assumed that 
all pions seen in the detector reflect the physics at 
the time of freezeout, but it is thought that roughly 
half of the detected pions come from the later decay 
of resonances |20]. Let us look at how this affects 
our results. Consider the peaks in Figs. [2] and [3j in 
the vicinity of fig where freeze-out occurs closest to 
the critical point. Because the cumulants (for ex- 
ample K47r) are extensive, our calculation of the nor- 
malized cumulants (for example uj^t^ — ha-k / {^-n)) 
would be correct if the experimentalists measur- 
ing (jJiT^ divide by the number of pions present at 
freezeout. By dividing instead by the number of pi- 
ons seen in the detector, the experimentalists will 
obtain a smaller ui^t^ than in our calculation. This 
is an effect that can be corrected for. 

• There are physical effects that can limit the up- 
ward fluctuation of Np. For example, if the pro- 
ton number density becomes too large, it will not 
be a good approximation to treat the protons at 
the time of chemical freezeout as noninteracting. 
In Appendix A we make a crude attempt to esti- 
mate the consequences of including such effects on 
the normalized cumulants. It will be much easier 
to model the consequences of this effect with data 
that show evidence for critical fluctuations in hand, 
since such data itself will indicate whether upward 
fluctuations in Np are cutoff, and if so at what Np. 

• The fact that net baryon number is conserved will 
also limit the fluctuation in Np. The magnitude 



of this effect depends on the size of the acceptance 
window, and for noncritical (and Gaussian) fluc- 
tuations has been studied in Refs. |31 It 
also depends on the features of baryon number fluc- 
tuations outside the acceptance window. It may 
translate into a sharp cutoff on the upward fluctu- 
ation of Np (as, e.g., proposed in a model study in 
Ref. |3^) or the reduction in flucutations may be 
more smoothly distributed over a range of Np. We 
defer investigation of this effect to future work. Ex- 
perimentalists will also be able to learn more about 
this and other effects by studying the dependence 
of the normalized proton cumulants on the width of 
the rapidity acceptance window, once there is data 
showing evidence of critical fluctuations. 

• We have focussed on fourth and lower order cumu- 
lants. Our results show, though, that higher order 
cumulants depend on even higher powers of the cor- 
relation length ^, making them even more sensitive 
to the proximity of the critical point. However, the 
measurement of higher order cumulants involve the 
subtraction of more and more terms, making it im- 
portant to determine the precision with which they 
can be measured. We have stopped at fourth order 
because current analyses show that these cumulants 
can be measured with small error bars. If cumu- 
lants beyond fourth order are measured, it will be 
possible to construct further ratios of observables 
that overconstrain parameters or are independent 
of parameters. 
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Appendix A: Toy Model Probability Distribution 



In Section II B we presented the calculation of (and re- 
sults for) the second, third and fourth cumulants of the 
pion, proton and net proton multiplicity distributions. 
We found, for example, that with our benchmark param- 
eters u!J^p peaked at a value of around 400, while for a 
Poisson distribution w^p = 1. This dramatic increase in 
LOip due to critical fluctuations with a reasonable value of 
C = Cmax = 2 fm raises the question of what the distribu- 
tion whose moments we calculated looks like — does it 
in any way look unreasonable or unphysical? Although 



the results from Section II B determine arbitrarily high 
cumulants of the proton multiplicity distribution, they 
do not allow us to determine the shape of the distribu- 
tion itself. In this Appendix, we provide an example 
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of a probability distribution ^^{Np) which has values of 
LUip for i = 2, 3 and 4 that are comparable to those we 
calculated in Section |II B[ This toy model distribution 
is somewhat, but not completely, ad hoc, since we shall 
construct it in a way that does reflect the origin of the 
critical contribution to the fluctuations. 

Let us consider a free gas of particles of a given species 
that have a mass M{a) which is a function of a back- 
ground field fj. As an example, for protons we shall use 
M((t) = rrip + Qpa. And, let us assume that the a field 
fluctuates with a given probability distribution V^icr). 
The central simplification that we are making in con- 
structing this toy model is that we are assuming that a 
is spatially homogeneous. The field a fluctuates, but at 
any given time it is the same everywhere in space. Let us 
model the probability distribution for the number Np of 
particles with mass M(cr) by considering the fluctuations 
of Np due to the fluctuations in a. Integrating over the 
fluctuating a we obtain the probability distribution for 




FIG. 6: An example of a distribution with ^ 400. 

The construction of the model distribution is described in the 
text, as are the values of its first few cumulants. Np is the 
number of protons in a volume V = {5 fm)^ in the toy model 
distribution. Other parameter choices are described in the 
text. 



V{Np)^ j daVa{a)PM(.){Np) (Al) 

where Pm{N) is the probability distribution for a particle 
with fixed mass M which we choose to be Poisson: 



Pm{Np) 



N^p 
AM 



-TV 



(A2) 



where N is the expectation (mean) value of Np for the 
distribution Pm(Np). In thermal and chemical equilib- 
rium, 



N = V 



1 



k exp ( ^^^!±MM!z^ I ± 1 



(A3) 



are 



where we choose the positive sign since protons 
fermions. N depends on M and, therefore, on a. 

The probability distribution of a is determined by the 
effective potential n{a): 



Va{(j) - exp -V 



T 



where the effective potential can be written as 



^ ' 2 



A, 



A, 



(A4) 



(A5) 



namely (2.2) without the spatial gradients. Eqs. (Al) 



(A5) define the probability distribution for the particle 



number, which will depend, among other things, on the 
correlation length ^ = m"^. Note that since the vol- 
ume V in the model corresponds to the volume within 
which the a field is homogenous we should think of V as 
a parameter in the toy model just as f is. The model 
treats only the zero-momentum mode a = ^^a{x)/V 
of the critical field, ignoring all other modes, i.e., the 



space variation of the field aix). This means, in par- 
ticular, that it ignores the fact that the correlations are 
exponentially small between regions of space separated 
by distances further than ^. For this reason we should 
not choose a value of T/^/'^ that is very much larger than 
t 

As an example, in Fig. [6] we plot the toy model prob- 
ability distribution for the number of protons, V{Np) of 
(Al), with ^ = 2 fm, A3 and A4 taking their benchmark 
values, hb = 400 MeV, V ^ {5 fm)^ and gp = 6.185. We 
can then evaluate the mean and cumulants of this toy 
model probability distribution, and we find 



□del 
□del 

2p 

model 
3p 

model 



5.2, 
4.5, 

37, 
405. 



(A6) 



We chose all parameters in the toy model at their bench- 
mark values with the exception of gp, whose benchmark 
value is 7. We chose gp = 6.185 in the toy model in 
order to get a probability distribution whose fourth cu- 
mulant is similar to that we calculated in Section IIIBI 
In our full calculation of Section II B[ with ^ = 2 fm, 
/is = 400 MeV, and all parameters at their benchmark 
values including in particular gp = 7 we obtain 



LL>2p 
^3p 
UJ4p 



3.0, 
4.2, 
30, 
405 



(A7) 



where we have quoted {Np) — Vup for y = (5 fm)'^. 
(The w's calculated in Section II are intensive, meaning 
that they are the same for any choice oiV .) We see that 
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the distribution in Fig. [6] has cumulants that are similar 
to those we calculated in Section II B[ including in par- 
ticular having as dramatically large a value of a;4p. We 
see from the figure that an uj^p that is ^ 400 times larger 
than the Poisson value does not indicate an unusual look- 
ing distribution. We can also note that the large positive 
4th cuniulant is a consequence of the skewness of the dis- 
tribution. If the distribution were symmetric, the large 
positive 4th cumulant would indicate a highly peaked 
distribution, but not so here. 

We can also use our toy model to make a crude es- 
timate of how our results for the cumulants would be 
affected by an upper cutoff on Np. In a heavy ion col- 
lision, Np (say in one unit of rapidity) cannot fluctuate 
to arbitrarily large values. We can implement this in the 
toy model by putting an upper cutoff on Np. Lets as- 
sume that the neutrons are fluctuating with the protons, 
as is in fact expected [30l. It seems clear that nucleon- 
nucleon repulsion (that we have not taken into consid- 
eration) would cut off upward fluctuations in Np some- 
where below those that correspond to nucleon densities 
of l/fm'^, meaning 60 protons per (5 fm)^ volume. To 
get a sense of the size of these effects, we tried cutting 
off the distribution in Fig. |6]at Np = 30. We find 



utoff 



□ del 



P/ 



cutoff model 
2p 

, .cutoff model 

, .cutoff model 



5.1 , 
4.4, 
33 , 
289 . 



We see that a cutoff like this has little effect on the 2nd 
and 3rd cumulants, but it does reduce W4p by 28%. (See 
also Ref. gD] for a study of the effects of introducing a 
cutoff at large Np in the absence of critical fluctuations.) 



Appendix B: Mean Transverse Momentum 
Fluctuations 



and similarly for k^^Spt)- We can normalize k^^Spt) by 
defining a dimensionless and intensive variable Fk'. 



{N)^-^ku{5pt) 



(B2) 



where (A^) is the total particle multiplicity and vf^^(pT) 
is the variance of the inclusive (single-particle) px- 
distribution, defined as 



Apt) 



1 



/k("k) 



(kT - prf ("k) 



(B3) 



Upon evaluating ( |B2[ ) using the correlators given in sec- 
tion ITlSI wc obtain the critical contribution to Fk- For 



pions with /i^ = at T = 120 MeV and £, = 2 fm, we 
find 



^3 



-0.0131 and F^^ = 0.0177. 



(B4) 



In addition to the critical point contribution, expression 
(B2) receives contributions from Poisson statistics, Bose- 



Einstein enhancement, resonances, effects of radial flow, 
etc. It was shown in [20] that the effects of resonances 
and radial flow are very small and hence we will ignore 
them. Here we compare the critical point contribution 
to that coming from Bose-Einstein enhancement. The 3- 
and 4-particle correlators for an ideal Bose gas are given 
by 



(AS) ((('5nk)^))BE = (?^k)((?^k) 



BE 



(«k)((n-k) 



l)(2(nk) + l), 
l)(l + 6(nk)((nk) 



(B5) 
1)),(B6) 



where here by (nk) we mean the mean occupation num- 
ber for an ideal Bose gas. Evaluating (B2) using the 



above correlators will give us the Bose-Einstein and the 
Poisson contribution to Fk- In order to isolate the Bose- 
Einstein effect we subtract the Poisson contribution (nk) 
from Kfc and then evaluate i^^. Using the same parame- 
ters as above we find 



The correlators found in section III Al can also be used 
to estimate the effect of the long wavelength fluctuations 
in the vicinity of the critical point on higher cumulants 
of the mean transverse momentum px- For example, the 
cubic cumulant of the mean px distribution around the 
all event mean pf, namely K3{5pt), is given by 

KsiSpr) = {{{PT-Prf}) 

= 7y~"TT3 / / / ([ki]T -PT)([k2]T -Pt) 

(jk("-k)) "'ki Jk2 Jk3 

X i[k3]T - PT){{SnkiSni,^672k^)) , (Bl) 



pBE 
-^3 



-0.2480 and Ff^ = 0.9388 



(B7) 



We see that the contribution of critical fluctuations is 
smaller than that due to Bose-Einstein effects. We con- 
clude that it would be very difficult to use higher cu- 
mulants of the mean px distribution in order to search 
for the critical point. Furthermore, as kinetic freeze-out 
(where particle momenta freeze) occurs after chemical 
freeze-out (where particle numbers freeze), it is easier for 
Pt fluctuations to get washed out (see, e.g. Ref. 41 ), 
making them even less favorable observables in searching 
for the critical point. 
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